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Abstract 

l/~, . We study a simple toy model, which although probably does not have any di- 

rect physical applications, can serve as a nice pedagogical example for explanation 
strange phenomenon of negative radiation pressure and can also give some insight 
for understanding the interaction between radiation and a vortex. The model dis- 
cussed is a classical field theory of two interacting scalar fields in 1+1 d. Without 
coupling one of the fields is governed by the ordinary <ji 4 equation and the second 
one obeys the Klein-Gordon equation. The possibility of existence of the negative 
radiation pressure with a respect to the mass ratio of those two fields is discussed. 



1 Intoduction 

>, 

In modern physics topological defects play a very special role. Because of huge variety of 
topological defects they found applications in many branches of physics starting from glueballs 
and fluxtubes in QCD through vortices in quantum liquids and superconductivity, ferromagnets 
to cosmological strings and domain walls. They owe their popularity to their stability due to 
f""^ ' nontrivial asymptotic behaviour. In many processes they can be treated as ordinary pointlike 

^^ \ particles (kinks, monopoles), strings or membranes. However they have also many properties 

when their internal complex structure becomes dominant showing that topological defects can 
be very different from particles or other simple objects. 

One of the phenomenon in which, on the first sight, topological defects behave like ordinary 
particles is radiation pressure. When a kink is hit with scalar radiation it reflects some of this 
scalar field. This reflection pushes the kink. This effect can be effectively described as radiation 
pressure very similar in its nature to radiation pressure known from electrodynamics. However 
there is a large group of topological defects which behave in a completely different manner 
revealing the negative radiation pressure (NRP) like kinks in </> 4 Q]|2] and vortices in Goldston's 
and abelian Higgs models. 

The NRP was first discovered in </> 4 model. When the plane wave hit the kink it was acceler- 
ating towards the source of radiation. This effect was confirmed using both numerical simulation 
and perturbation theory. In the first order the kink is transparent to the wave due to the so called 
reflectionless potential. Therefore our usual intuition can be missleading. The further corrections 
showed that after the transition the wave carried the same amount of energy but more momen- 
tum. The surplus of momentum forces the kink to accelerate towards the source of radiation. 
The force is proportional to A 4 . In other models which do not have kinks with reflectionless 
potentials the solitons behave as one could expect. They are pushed by radiation and the leading 
term in the force is proportional to A 2 . 



Here we propose a different example of the negative radiation pressure. We have decided to 
discuss a strange model because of a few reasons. The appearance of the NRP in this model has a 
different mechanism. The leading term in the force is proportional to A 2 and therefore it suggests 
that it comes in the first order of perturbation scheme. Moreover the system possesses another 
parameter (a mass of Klein-Gordon field) which one can change and find two different regimes: a 
negative radiation pressure and the positive one. Our preliminary numerical calculations suggest 
that similar mechanism is responsible for the NRP in case of vortices. The NRP is proportional 
again to A 2 . Moreover in the abelian Higgs model there is also a parameter which can be 
triggered. It is usually used to distinguish between superconductors type I and type II. This 
parameter possesses also a different value which distinguishes between vortices exhibiting NRP 
and PRP 

The model discused in the present paper has very similar properties of the radiation pressure 
as Goldstone's model. The adventage of the model discussed is its simplicity. 

The paper is constructed in the following way. First we give a motivation as to why we have 
chosen this particular model. Next section is a description of the model itself and some spacial 
solutions existing in this model. In the following section we analyse small perturbations around 
topological defect using standard perturbation procedure both in coupling constant and wave 
amplitude. Within this perturbation scheme we construct a wave travelling from rhs of the kink. 
We derive an analytic form for the force exerted on the kink by the wave. The next section is 
devoted to numerical results where we compare our perturbation scheme results with numerical 
solution to the full partial differential equation. Finally we summarise the results and discuss 
possible applications and extensions. 

In the Appendices we give some more detailed calculations. The last Appendix is a short 
description of Goldstone's and Higgs modeles pointing out the similarities to the model presented 
in this paper. 

2 Motivation 

As we mentioned before the negative radiation pressure was first observed numerically and then 
explained analytically in case of one-dimensional kinks in (j> 4 model. The key to this unusual 
phenomenon is transparency of the kink in linear regime. The second order calculation showed 
that a wave with twice the frequency is created and that wave carries more momentum than 
the incident wave creating the surplus of momentum behind the kink. Because of the energy 
and momentum conservation laws the surplus of momentum must be compensate with kink's 
motion. The kink starts moving towards the source of radiation. This is purely nonlinear effect. 

The feature of reflectionless in linear regime is rather rare and does not necessarily lead to the 
negative radiation pressure. As a counter-example in paper (2) we have discussed another often 
studied model i.e. sine-Gordon. Kinks in this model are also transparent but they are transparent 
exactly to all orders in perturbation series and hence experience no radiation pressure at all. 

Fortunately the phenomenon of NRP turned out to be quite robust with respect to pertur- 
bation of the (f: 4 potential. Although the reflectionless was lost the nonlinear effects were still 
visible and above a certain critical amplitude the NRP was present even for quite large deforma- 
tion of the potential. 

As we stated in the introduction the NRP is visible not only for kinks in 1+1 d but also 
for Goldstone's and abelian Higgs model vortices. Unfortunately the analysis of the negative 
radiation pressure is difficult for those models both analytically and numerically and will be 
presented in separate papers. We believe that in case of vortices the negative radiation pressure 
has a different nature than in <f> 4 model. Our preliminary numerical analysis shows that in case 
of vortices the force exerted by travelling wave is proportional to A 2 on the contrary to A 4 in 
<t> 4 model. Moreover the vortices are clearly not reflectionless and hence a different mechanism 
must take place. 

In the present paper we propose a model exhibiting kinks undergoing negative radiation 
pressure proportional to A . We explain the mechanism which stands behind this process which 



is different than the one in 4 model. The model is quite simple and probably does not have any 
direct physical application but some of the features seem to be generic and are not blurred with 
technical difficulties. 

When the negative radiation pressure is present there always must be some kind of mech- 
anism creating a surplus of momentum behind a defect. In </> 4 case the momentum to energy 
ratio was k(cj)/uj = \/cj 2 — 4/oj. The wave with twice the frequency has momentum to energy 
ratio equal to k(2aj)/2io = \/u! 2 — 1/w > k(u)/u). Therefore a wave with twice the frequency 
created due to the nonlinearities carries more momentum than the incoming wave. 

This is not the only possible way to create a surplus of momentum. Let us consider two 
scalar fields with masses mi and 7712 scattered on some object (i.e. topological defect). Suppose 
that only a wave with mass mi comes from one side of the defect. During scattering process 
the wave is converted into a wave with mass 7712. The momentum to energy ratio is equal to 
y/ui 2 — rrii/uj and therefore if 7711 > 777,2 we have again a surplus of momentum behind the 
object which could undergo the negative radiation pressure. Similar process is possible for a 
single scalar field with degenerate vacua when the excitations around those vacua have different 
masses. 



3 Model 

Let us now propose a simple model of two interacting fields with two different masses. We want 
the model to exhibit topological defects. A natural candidate would be a 4> 4 model field coupled 
with some other fielqj: 



C = i (</>? - &) + \ [f t - &) - i(V - 1) + ^ - \{m 2 


~K 2 )i? (1) 


or equivalently 




£ = 2 ( 0? - 0* - O 2 - I) 2 J + 2 (^t ~^l- m 2 ip 2 ) - k(4> 2 - 


- l)tf, (2) 


where m is a mass of tp field and k is a coupling constant. Note that when « 2 = 


; m 2 the vacuum 



vacuum consists of merely two disconnected points [cj> — ±l,tp = 0) revealing Z2 symmetry 
and domain walls (in general) or kinks (in one-dimensional models) are possible. 
The equations of motion are: 

<t)tt - (f>xx + 2<i>(</> 2 - 1) + 2k0V = 0, 
i>tt — 4>xx + m 2 i/j + k(4> 2 — 1) =0. 

Unfortunately we do not know even nontrivial time independent solution - the kinks in analytical 
form. Numerically calculated profiles are depicted in FigureQ] For small values of k the cf> field 
is very similar to tanha; as in ordinary (f> 4 model. The ip field is a gaussian-shape function 
with its height proportional approximetly to «. Large x value asymptotic approximation gives 
4> ~ 1 — ae~" x and ip ~ e~" x with n being the smaller value out of: 



+ 4±^/(m 2 -4) 2 + 16K 2 



'Other possible models are discussed in section |5~T1 



Kink profile for re =0.1 and m— 3.0 
Fieldd 




Figure 1 : Static solutions for the kink for k = 0.1 and m = 3.0 



4 Radiation pressure 
4.1 Planar wave 

Let us consider a process when a kink is hit with a ?/>-field wave with frequency u) coming 
from +oo. In order to calculate the leading term to the radiation pressure we do not need to 
know the full kink solution. All we need is the the perturbation series in coupling constant k: 

4> = 4> {0) + k<P w + k 2 (2) + • • • and V> = ip {0) + K.ip {1) + kV 2) + • • • • 
In O(k ) the equations take the form: 



*£? 



■1>, 



if,™ - Vi°J + m> 



-lj=0. 

2.;.(0) _ q 



(5a) 
(5b) 



For static kink solution from these equation we obtain 4>^°' = tanhx and ip = 0. In 
order to study a scattering process instead of ip = we take ip —A cos(kx + ujt) where 
A is an amplitude of the wave and k — \/lo 2 — m 2 is a wave number. The first order equations 



4l } - £1 + 2 (3^ 2 - l) 0« + 2*< V 0) = 0. 



-1=0. 



(6a) 

(6b) 

The equation J6bt has only static inhomogeneous part so it gives only correction to the static 
solution of the kink. Since we know solutions to the homogeneous part (for ipi — 0) we can 
express the solution of inhomogeneous equation in terms of Green's function: 



i (i) ^ 



dx 



+ e 



: / dx' 

J X 



cosh 2 x' J x cosh 2 x' 

These integrands can be calculated using hypergeometric functions. 



(7) 



Since ip(°> = l e lu}t + kx -\-c.c. the solution to the equation d6at can be sought in the fallowing 

form: 

(i) = l e i^ +(i)+ l e -i^ (a;) (g) 

where £+ is a solution to 

\ dx 2 cosh 2 x 



f—+q 2 + l 2 - )$+ = 2Ae ifc *tanhx, 



with q 2 — uj 2 — 4 and £_ = £+. Fortunately we already know the solution r\ q to the homoge- 
neous part of the above equation (see l2ln 

, . 3 tanh 2 x — 1 — q 2 — 3io tanh x iax ,,_, 

r)q(x) — ; e^ (10) 

N/(« 2 + l)(? 2 + 4) 

and again we can use the Green function technique to construct the solution: 

f+(x) = -y- f dx' n q {x')f{x') - ^ |V r?- 9 (x')/(x'), (id 

where W 7 = — 2ig is Wronskian and f(x) = J 4e li:x tanhx is the rhs of the equation J9). As for 
all scattering processes we are interested only in the asymptotic form of this solution. Therefore 
it is convenient to rewrite the above solution for large values of a; in a little different form: 

rn poo -i roo 

£ = --=2-1 dx' n q {x')f{x') + — J dx f(x') ( V - q (x)T] q (x') - Vq (x)rj- q (x')) . (12) 
To calculate the second integral we can use the asymptotic forms both of r) q (x) 

ih (a ->(»)» 2 -g 2 - 3i 9 =e ia» (13) 

^(q 2 + l)(q 2 +4) 
and /(x) = 2Ae ikx to obtain: 

I = — / dx' /(a; ) (»?-a(x)?j,(x') - %(ar)»/-,(x')) = 



— / dx' I exp(i(gx' — qx + kx')) — exp((i(— qx' + qx + kx')) 
"* J x V 



(14) 



With help of identity 



dx e ax = nS(a) H — e ax (15) 



a 
we obtain 

I = 2 k^ + ^w lm+q) - 5{k - q)) - (16) 

Because only when m — 2 (when model does not exhibit topological defects) the <5s give contri- 
bution, I represents only the inhomogeneous part of the solution. Note that the same solution is 
present for x — > — oo. 

The first integral of equation d 12b can be calculated via residua method. The solution for 
large values of x can therefore be written as: 

ti + (x^+oo)=R{q,k)r,- q (x)- J_ - . (17) 

where 

, v T (3k 2 - q 2 - 4) 

R(Q, k ) = V = ' , tt : (18) 

2«V(9 2 + 1 (« 2 + 4 sinh (2±^tt) 



2 Note that there is no reflection part proportional to e ' 9 



is reflection coefficient. For x — > — oo similar procedure leads to the following solution 

e+(a; -► -oo) = r(g, k)r,-q(x) + j£— 5 . (19) 

where T(g, fe) = R(—q, k) is the transition coefficient. One can see that for q = fc <4> m = 2 
the transition coefficient becomes infinite. Of course this means that our perturbation scheme 
fails in the vicinity of this point. 

Another important observation is that for q 2 = 3k 2 — 4 <4> uj 2 — ^j- the nominator vanishes 
and the kink becomes again. This should be clearly seen in numerical simulations. 

4.2 The force 

Having the asymptotic form of solution representing travelling wave 

il^we 1 "' I R(q, k)r]- q (x) — fc g_ 5 1 + c.c. for x > +00 
■ } lkx\ ( 20a ) 

|j4fi;e lT(q,k)r/ q (x) + ,|_ a 1 + c.c. for a; > — cxj 

I ^4cos(wt + kx) for 2 — > +00 
ip(x,t) = { \ \ (20b) 

I Abcos(u>t + kx) for x — > — c» 

we can calculate the force which is exerted on the kink by this wave. In the second equation we 
changed the amplitude from A to Ah. Actually as we will shortly show b = 1 + 0(n 2 ) so this 
correction is not visible at this point of perturbation series but it is necessary to fulfil the energy 
conservation law. If the kink is initially not moving the rate of energy flowing into a large box 
containing the kink is given by (for more details see Appendix lAl 

{E) T ={<l>t<t>*+i>t'>Px)T\ , (21) 

were we average over a period. After substitution we obtain: 

^\ =^LoA 2 (-qK 2 R 2 + k-kb 2 -qK 2 T 2 ) (22) 



and since energy must be conserved ((E) — 0) we obtain 

'■1 
k 



K 2< t(R 2 + T 2 ) . (23) 



note that our scheme is obviously invalid when b < i.e. 

R 2 + T 2 >-^-. (24) 

K^q 

This correction to the amplitude of ip wave would be obtained after calculating next step in our 
perturbation series but since the energy must be conserved we obtained this result as a consis- 
tency condition. (Similar result was discused in more details in (2).) 

The force exerted on the kink can be calculated using the momentum conservation law(for 
more details see Appendix [At: 

F =/^\ = ^A 2 (-k 2 - K 2 R 2 q 2 + b 2 k 2 + k 2 T 2 ) (25) 

or using l l23l : 

F=±A 2 q(T 2 (q-k)-R 2 (q + k)). (26) 

The above force is proportional to A 2 (contrary to A 4 in (f> 4 model (2)). Moreover it can both 
be positive (kink accelerates towards the source of radiation - negative radiation pressure) and 
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Figure 2: Force exerted on a kink (divided by k 2 A 2 ) for three different values of m. 



negative (the kink accelerates with the wave due to the positive radiation pressure). Note that 
when q < k i.e. m < 2 the force is always negative (or at most zero) whatever the coefficients 
R and T are. When m > 2 the direction in which the kink accelerate can be determined only 
after substitution the values of R and T. The negative radiation pressure appears when 



T 2 > g±k R 2 

q — k 

This inequality can be rewritten (when R 7^ 0) as 



sinh 



g— fc. 



sinh 



9 + fe, 



> 



(27) 



(28) 



which is true for all q > fc > 0. The force vanishes only when R — that is when a; 2 = '^—. 

Another important conclusion follows from our considerations which distinguishes the me- 
chanizm presented in this paper from the NRP in (j) 4 . The force in our model comes directly 
from the first order calculations in amplitude of the wave. We do not need therefore to limit our 
consideration to monochromatic waves. Any i/j-wave can be expressed as a wave packet: 



/ 



ip(x,i) = /dfe A(k)e i 



(kx+ui(k)t) 



+ C.C. 



The averaged force exerted by the packet is therefore: 



T = -k jdk A 2 {k) (T 2 {k)(q(k) - k) - R 2 (k)(q(k) + k)) 



(29) 



(30) 



Because when m > 2 the force is not negative for all frequencies and the kink would always 
undergo the NRP. In case when m < 2 the kink would be pushed by the packet. 



5 Numerical calculations 

As we have stated before, even in this simple model we do not now the exact analytical form 
of static solutions. Therefore they need to be obtained numerically. We have done this using 
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Figure3: Trajectory of the kink along with fitted parabola X(t) = \at 2 +vt + c for k = 0.1, m = 2.5, 
A = 0.1 and w = 3.5. The fitted values: a = (7.57898 ± 0.00069) ■ 10" 6 , b = (3.23 ± 0.18) ■ 10~ 5 
andc= (1.1 ±1.9) ■ 10" 4 . 



collocation Chebyshev spectral method in variable s = tanha;. In this way we obtained the 
solutions 4>o(x) and tpo(x) depicted in the FigureQ] 

The next step was to find the initial condition representing travelling ip wave. Obtaining a 
solution satisfying all assumptions of the travelling wave is rather cumbersome. Therefore we 
have decided to use only an approximation: 



(x,t — 0) = (/>o(x), 



,t = 0) = 0, 



(31) 



and 



ij)(x,t — 0) = ipo(x) + Acos(kx), ip(x,t = 0) = — Aujsm(kx), (32) 

Disadvantage of the above initial data is that at the very beginning the kink has a no-zero ini- 
tial velocity. This is due to the fact that the above initial conditions are combination of many 
eigenstates of the linearised equation around the kink, including translational mode. However 
this type of conditions show how roboust and generic our considerations are. 

We solved our PDE equation (|3j using five-point descritization of the second spatial deriva- 
tive and integrating the obtained system of ODEs with 4th order Runge-Kutta method. We have 
approximated the position of the kink with the topological zero of the cj> filed. The position of 
a static kink is well defined and it coincides with the topological zero. In case when the kink is 
exposed to the wave its position is not well defined but it cannot be far away from topological 
zero. During its motion the kink recedes from its initial position significantly and those subtleties 
can be neglected. The example trajectory of the kink is presented in the Figure [3] The path of 
the kink is very similar to parabola (as expected) and we can fit it with a function in the form 
X(t) = | at 2 + vt + c obtaining its initial velocity v and acceleration a. 

In the Figures [4] and [5] we have depicted the measured acceleration as a function of mass 
parameter m and frequency u along with theoretically predicted functions (force from equation 
J26b divided by kink mass for k — 0: Mk = f )• Figure [4] shows that analytical predictions 
are quite precise (outside the shaded region - vicinity of m = 2, when our perturbation scheme 
fails). Note that the force disappears (kink is transparent to the radiation) when 2u> 2 — 3m 2 . 
This can be seen in both Figures. Our analytical prediction differs significantly from numerically 
obtained acceleration also for small frequencies (Figure [5). This can be due to the fact that 
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our approximated initial conditions differ significantly from the travelling (eigenfunction) wave, 
especially for small frequencies. 

Except this low frequency limit and vicinity of m = 2 our theory conicides with numeri- 
cal simulation with average error of a few percents which is satisfying having in mind all the 
assumptions and simplifications we used. 

The last Figure shows how acceleration grows with amplitude A and a coupling constant 
n. From our theoretical prediction we expect that acceleration should be proportional to A 2 k 2 . 
In the Figure[6]we have plotted the acceleration as a function of An. There are two lines: one 
for constant A — 0.1 and the second one for constant k = 0.1. Those functions coincide very 
well especially when both k and A are small. We have fitted the exponent for both functions in 
the form a = ae^ K . For constant k we obtained /3a = 2.0888 ± 0.0075 and for constant A: 
(3 K — 2.0196 ± 0.0025. This indicates that indeed the acceleration is proportional to A 2 k 2 . 



5.1 Other models 

5.1.1 Other two-component fields 

We have addressed the question of the negative radiation pressure also for other models. Namely 
we have tested models of similar structure as in equation {2j: 



Ul - 4>l - (0 2 - l) 2 ) + \ (ft -ill- m 2 ^ 2 ) - Vi„t(<j>,TP), (33) 



£ = |(« 



where Vint (0, ip) is an interaction potential between fields </> and ip such that the interaction is 
weak far away from the kink, and large very close to the topological zero. In the present paper 
we have focused on potential: 



Another potential we have used had the following form: 



(34) 



v£>(<t>,i>) = K(4> a -iyil>. 



(35) 
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Although this potential defines theory with unbounded energy form below, the system is stable 
with respect to enough perturbation around the kink solution. To fix that inconvenience we could 
add some other potential which is negligible for small values of field and dominating for large 
values. 

Other potentials (giving energy bounded from below) were: 

V&W) = *fa a -l)^ a . (36) 

^WH^-OV. 07) 

For values k = 0.2, A = 0.2, u) = 2.5 and m = 2.2 we have performed numerical simulation 
and we have found that apart from V> n \ only for potential V!> n \ the kink revealed the NRP. For 
the other two potentials the kinks were pushed by the radiation. 

5.1.2 A halfcompacton 

We have observed an interesting behaviour of a kink in a model of one interacting field with 
asymmetric potential described by lagrangian density: 

C=\{ti-$l-\^-m+lf). (38) 

In this model the potential has two equal minima <j) — ±1, therefore topological defects are pos- 
sible. Important difference from other well known models is that the filed of the kink approaches 
vacuum (</> = 1) at finite x while on the other side it approaches vacuum exponentially. This 
type of object is usually referred to as a halfcompacton. Compactons are topological defects 
which differ from vacuum only in a finite segment [6J. They appear usually when the potential 
is V-shaped around vacua. 

In this particular model the potential is smooth around one vacuum {<j> = —1) and sharp 
around the second one. Small perturbations around <j> — — 1 can be described as field with mass 
m 2 = 2. Formally, the second derivative of the potential at the second vacuum is infinite and 
so is (formally) mass of small perturbation. Arbitrary small perturbations around the V-shaped 
potentials are nonlinear. 

We have found that the kink is pushed when is hit with radiation coming from the vacuum 
<t> = — 1. However it reveals NRP when hit with the wave is coming from the second vacuum 
<jf> = l. The halfcompacton is always pushed towards the V-shaped vacuum. Again this process 
can be understood because of the mass difference {m\ = 2 and formally m?, — oo) of small 
excitation around two different vacua. 



6 Conclusions and discussion 

In the paper we have presented a model of two interacting scalar fields. Without interaction one 
of the fields (cj>) obeys cj> 4 and the second (i/>) Klein-Gordon equation with mass m. This model 
reveals Z2 symmetry and topological defects, like kinks, are possible static solutions. We have 
shown that the kink, when exposed to scalar radiation of field ip, can undergo a strange phe- 
nomenon of the negative radiation pressure. This happens when the mass parameter m is larger 
than the mass of <j> field which in our scaling is equal to 2. Below that limit the kink is a subject 
to positive radiation pressure. We have calculated analytically the force exerted on the kink by 
the radiation in the leading order both in amplitude of the wave and coupling constant. The force 
was proportional to product of squared amplitude and coupling constant. The results were con- 
firmed with numerical simulation of the full partial differential equation. The difference between 
analytical predictions and numerics were caused by simplifications used during calculations and 
preparation of initial conditions for numerical simulations. The descrapency was mostly a few 
percent. 

The model can serve as an example of a different mechanism of negative radiation pressure 
than the one presented in our previous papers (T|[2l in the case of <f> 4 model. The NRP in (f> 4 is 
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a purely nonlinear effect and the force is proportional to the the fourth power of amplitude. The 
case of (p 4 is a special one, due to the transparency of the kink, which is a rather rare feature. 
In other models kinks undergo usually positive radiation pressure. However, it is possible that 
there can be a critical amplitude above which the nonlinear effects, like NRP, are dominating. 
Another problem with the 4 is that it is difficult to say how the kink would react to the non- 
monochromatic waves. As we have stated before the effect is nonlinear and no superposition is 
possible. 

The NRP presented in the present paper has different properties. It emerges in linear approxi- 
mation and therefore a superposition rule can be applied to it. This mechanism seems to be more 
generic than in the <f> 4 model. We did not need the transparency of the kink (in linear regime). 
The negative radiation pressure, although a very strange phenomenon, is probably more common 
than we had initially expected. We have changed the interaction potential between fields <j> and 
ip. Then we performed numerical simulation. Some models revealed only positive radiation 
pressure. However, we have found a few models with NRP. For all of them there was a critical 
value of mass parameter m cr = 2 above which we observed NRP and below PRP The observed 
acceleration was proportional to A 2 . From whole that class the model presented in this paper 
was the simplest one. The simplicity allowed us to perform analytical calculations, which could 
be difficult in other models. 

Our preliminary numerical simulations show that also vortices in Goldstone's and abelian 
Higgs model reveal NRP which is in its nature very similar to the one presented in this paper 
(see Appendixlct. It is also possible that it appears in other more complicated models, especially 
when there are many interacting fields. Possible candidates can be for instance models of kinks 
studied in J4|- 

Importance of the NRP is still an open question. However, it may be important in many 
physical processes where topological defects are present. For instance it may be important for 
stability of a system of topological defects. Any small excitation would radiate being a source 
of long-distance attractive interaction. This issue is clearly dimmension dependent. In 1+1 d the 
defects accelerating towards the source of radiation would eventually collide and bounce back 
or annihilate usually creating a long-living metastable object called oscillon J5J - These oscillons 
can also play important role in collisions with other objects. In higher dimensions the collisions 
are usually less frequent and the radiation pressure may be more important. 

A Noether's theorem 

Our model possesses both time and spatial translational symmetry therefore we can define ener- 
gy-momentum tensor using the standard form: 

dd^cp ddpip 

where as usual T 00 is energy density and — T 10 = — T 01 is momentum density respectively. 
The elements of this tensor are: 



T°" =£ = \\< 



U 2 + & + (0 2 - l) 2 ) + i (i/> 2 + i? x + m 2 ^ 2 ) + K(<f> 2 - 1)V (40a) 
-V = <j>t<f>x + ipti/ix (40b) 



and 



T 11 =-U<8+4>l- {4? - l) 2 ) - \ (Vt +^l~ m 2 i> 2 ) + k^ 2 - l)t/>. (40c) 

Let us put a kink into a large segment [-L, L]. Using the Noether theorem we can find the 
rate of energy and momentum density flowing into this segment using only asymptotic forms of 
travelling waves. The Noether theorem states that 

d^T"" = 0. (41) 
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In particular, we can write 

d t £ = d T°° = -d x T w (42) 

which after integration over the whole segment gives the energy carried by the wave into the 
segment 



dE_f 

dt ~ J_ 



= / da; d x T w = frfa + M, 



(43) 



Moreover we also changed the sign to obtain the energy which flows into the segment and not is 
carried away by the wave. In the very similar way we obtain the ratio of momentum flowing into 
the segment: 



AP f 



= dxd T LL = 



-\Ul + <t>l- (4> 2 - 1) 2 ) - \ {jlh +i>l - mV) + <^ 2 - i)i> 



(44) 



Since we study interaction between a kink and a travelling waves, far away from the kink the 
filed has the following form (compare with equations d20all20bt ): 

4>(x, t) = ±1 + c\ cos(ujt + qx + Si) + C2 cos(ujt + kx + 62), (45a) 



if>(x, t) = C3 cos(ujt + kx). 
Substitution of the above into d43b and J44l > and averaging over time gives: 

L 



\dt) T 



-Ldqc\ + -ukcl + -ujkci 



(45b) 



(46a) 



As one can see from d20at the amplitude of second term in <j>(x, t) is the same on both sides of 
the kink and therefore they cancell each other out. In the similar way we obtain 



/dP 

\ dt 



I22 1,22 1,22. 

~2 q Cl ~ 2 3 ~ 2 ° 2 KC2C3 



(46b) 



Note that form equations d20all20bb and following discusion we can see that terms with C2 either 
have the same values at both ends of a segment or they differ by at least of order 0(k 4 ) therefore 
they do not give contribution to equations l !22t and l !26l l. 



B Asymptotic form of linear solution 

In our perturbation scheme we assumed that the coupling constant is small. We have obtained 
the solutions which asymptotically are waves with wave number equal to k = \fuj 2 — m 2 and 



^ 



4. It is easy to show that in the linear regime, far away from the kink, the wave 



numbers are different. This difference should appear at some point further in our perturbation 
scheme. Let us seek for solutions in the form 

4>{x,t) = cj> (x) + Ae^tix), 4, = V>o(a;) + Ae ilJt X (x) (47) 

assuming that A is small. After substituting the above into the equations of motion we obtain: 

d 2 



dx 2 



+ 2(30o - 1) + 2kVo 
2kcj>o 



'Ikx. 



v 



dx 2 













'£" 


2 

= LU 


~z 




A_ 




\_ 











(48) 



Far away from the kink, for x — > ±00 ipa — > and . 
form: 



±1. In this limit the equations take the 
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The eigensolutions are: 



dx 2 
±2k 



+ 4 



±2k 

_d 2 

dx 2 



2k, 

l 4-Lu 2 + k 2 . 2 

1 



where 



; 2 
«1,2 = 



+ 4 - 2cj 2 ± \/(m 2 - 4) 2 + 16k 2 



For small values of k 



and 



■ 2 2 2 . K 

k\ — lu — m — 4 — — 



fc 2 2 =a, 2 -4 + 4- ' 



+ 0(k 4 ) 



+ 0(* 4 ). 



(49) 

(50) 

(51) 
(52) 
(53) 



wr — 4 

The above results could be used as a test of how accurate our expansion in coupling constant is. 
One can see that the first correction to the wave numbers is of order 0(k 2 ). Note that for m not 
much different from 2 our expansion fails. 



C Application to vortices 



In this appendix we shortly discuss application of the mechanism described in the present paper 
to more physical, yet more complicated models. 

C.l Goldstone's model 

Let us consider a complex scalar field described by the lagrangian: 



c = d^&tf - (cj>*(f> - iy . 

This model reveals global U(l) symmetry. Vacuum manifold is a circle \<f>\ 
trivial solutions with cylindrical symmetry (vortices) are therefore possible, 
conciderations to cylindrical symmetry. The vortex solution has a form 

</>M) = /(r)e iArs , 



(54) 

= 1. Static non- 
Let us limit our 



(55) 



where N is a winding number (in this model only vortices with N — 1 are energetically stable). 
Perturbations, both around the vacuum state and a vortex, are of two types. So called Goldstone's 
mode is an excitation along the vacuum manifold (phase changing) and has mass equal to 
(m g — 0). The second perturbation excites the field perpendicularly to the vacuum manifold 
(changes field's amplitude, therefore we will call this excitation an amplitude mode). This mode 
has mass equal to 2 (m a = 2). 

Our preliminary numerical simulations show that vortices are pushed by planar wave gen- 
erated from Goldstone's mode. Amplitude wave exerts negative radiation pressure. The motion 
of the vortex is more complicated because vortices in this model are not well localized and have 
divergent energy. The phenomenon of the negative radiation pressure in case of vortices can be 
understood qualitatively using the same arguments as in the present paper. Because of the mass 
difference the amplitude mode creates a surplus of momentum behind the vortex. Numerical 
results along with analytical calculations in this case will be presented in a separate paper (3). 
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C.2 Abelian Higgs model 

Another very often studied model is abelian Higgs model for which the lagrangian has a form: 
C = D,4>(D^Y - ^ {4>*4> - rf) 2 - \F, V F» V , (56) 

where D M = <9 M — leA^ and_F M „ = 9 M A„ — <9„yl M is a Faraday tensor. Vortices are possible static 
solution due to the (local) C/(l) symmetry. Vortices in this model are more localized than those in 
Goldstone's mode. Again, our preliminary numerical simmulation show, that they are a subject 
to the NRP when hit with amplitude wave. Moreover the NRP also appears when the vortex 
is exposed to vector-field wave with appropriate polarization. We belive that this phenomenon 
has exactly the same mechanism as the one presented in the present paper: the mass difference 
between those fields (amplitude mode, phase mode and vector field) is responsible for this effect. 
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